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We present a compelling response of a low-dimensional strongly correlated system to an external 
perturbation. Using the time-dependent Lanczos method we investigate a nonequilibrium evolution 
of the half-filled one-dimensional extended Hubbard model, driven by a transient laser pulse. When 
the system is close to the phase boundary, by tuning the laser frequency and strength, a sustainable 
charge order enhancement is found that is absent in the Mott insulating phase. We analyze the 
conditions and investigate possible mechanisms of emerging charge order enhancement. Feasible 
experimental realizations are proposed. 

PACS numbers: 71.10.Fd, 74.40.Gh, 78.47.J-, 78.20.Bh 
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Nonequilibrium processes in strongly correlated elec- 
tron systems can provide new insights into the dynami- 
cal properties of these systems, which, in many aspects, 
can be qualitatively different from their weakly interact- 
ing counterparts. One such example is nonequilibrium 
induced phase transition [3-Q. As the system is driven 
away from the equilibrium, under certain conditions, a 
"crossover" from one state to another (metastable) state 
may occur. 

A well-known example is the insulator-to-metal tran- 
sition induced either by strong electric field or transient 
laser pulse, as a result of photodoping ■ In the case of 
Mott insulators, photocarriers are doublons and holons. 
In one dimension, just after the doping, generally, the 
system is in a metallic state due to the existence of itin- 
erant carriers and the benefit of the spin-charge separa- 
tion [H, [|[ . Experimentally, the dielectric breakdown in 
one-dimensional (ID) Mott insulators has been observed 
long ago in organic materials fiol . [Tl| and oxides fl2| . 
Photoinduced metallic state by transient light has also 
been investigated d, ■ More recently, the nonequilib- 
rium and nonlinear phenomena have stimulated signifi- 
cant interest within the cold atom community, where a 
novel realization of the fermionic Mott insulating state 
has been achieved flil ]. 

In this Letter we address the question of whether it is 
possible to reach other characteristically different states 
besides the metallic one after the photodoping. One 
possibility is the enhancement of charge-order when the 
attractive interaction between doublons and holons is 
incorporated [l5|. We seek such enhancement in the 
extended half-filled ID Hubbard model with an addi- 
tional nearest-neighbor interaction. At half filling, its 
phase is well-known and understood [T6l-[l9j. In large 
on-site Coulomb interaction, the model possesses two 
phases: spin-density-wave (SDW) and charge-density- 
wave (CDW), connected by a first order quantum phase 
transition, with algebraic decay of spin correlations and 



long-range charge order, respectively. In general, as a 
correlated system approaches the phase boundary, the 
response to an external perturbation becomes more elab- 
orate. In this Letter we show that in the latter case a 
substantial change in the electronic structure can be trig- 
gered by the optical pulse. In particular, when the system 
is originally in the Mott-insulating phase but close to the 
transition to CDW, a sustainable enhancement of charge 
order parameter is achieved by exposing the system to 
an external optical pulse with a rather carefully tuned 
frequency and amplitude. 

We consider the ID extended Hubbard model at half 
filling. The laser pump is incorporated by means of the 
Peierls substitution in the Hamiltonian: 



H(t) = -t h J2U A{t) cl a c t+ha +R.c. 



(i) 



where c| CT (c, ier ) creates (annihilates) electrons with spin 

, rii = 7ii,t + n i,h t h is tne hop- 



a at site i, rii 



ping constant while U and V arc the on-site and nearest 
neighbor interaction strengths, respectively. We model 
the external laser pulse in the temporal gauge via the 
time-dependent vector potential A(t) [2(| 

A{t) = A e-t t ~ t °W 2 % cos [co pump (t - t )] , (2) 

where Aq controls the laser amplitude, which reaches its 
full strength at t = to; td characterizes the duration time 
of light action. Notice that due to finite td, the incom- 
ing photon frequency is broadened into a Gaussian-like 
distribution, with the variance of 1/f ^ around the central 
value Wpump- We set t^ and t^ 1 as energy and time units. 



Wc now give a short overview of the equilibrium prop- 
erties (with A = 0) of the model, related to our work. We 
set U = 10 and vary V. With increasing V the system un- 
dergoes around £7 rs 2V a first-order phase transition di- 
rectly between SDW and CDW phases. Choosing rather 
large £7 we avoid the intermediate bond-order- wave phase 
that separates SDW and CDW at smaller £7 (l7Hl9j. In 
both phases charge excitations are gapped. In the Mott 
insulating, i.e., SDW phase, gapless spin excitations exist 
and the system displays no charge order. In contrast, a 
charge-density wave is characteristic for the CDW phase 
while spin excitations are gapped. 

In order to solve the time-dependent Hamilto- 
nian H(t), starting from the Schrodinger equation 
idip(t)/dt = H(t)ip(t), we employ the time-dependent 
Lanczos method, which is originally described in Ref. 1211 
and later applied and analyzed in more detail in Ref. [221 ] , 
followed by its applications in nonequilibrium dynamics 
of strongly correlated systems in Ref. [23[ . The basic idea 
is that we approximate the time evolution of \ip(t)) by a 
step- vise change of time t in small increments St. At each 
step, the Lanczos basis with dimension M is generated 
resulting in the time evolution 



\4>(t + St)) ~ e- iH V st \i)(t)) 



M 

i=i 



-ieiSt 



. (3) 

where e/ and | </>;), respectively, are eigenvalues and eigen- 
vectors of the tridiagonal matrix generated in Lanczos it- 
eration. We have confirmed that in our calculation, gen- 
erally M = 30 provides adequate accuracy when St < 0.1. 
For even smaller St, M can be smaller. 

In the succeeding numerical calculations, we employ 
periodic boundary conditions and set td = 5 and to = 
12.5. To investigate the time evolution of the charge 
order, we define the charge-charge correlation as 



C(j;t) 



(-iy 



L-l 



■i+3 



i)(m-i)m)), (4) 



where L is the lattice size. 

Figure [T] shows the results of charge-charge correla- 
tion functions C(j; t) for the system with 14 lattice sites, 
where the largest distance between two sites is 7. In or- 
der that the photon energy can be effectively transmitted 
to the electrons, we match the pumping frequency w pump 
with the first optical absorption peak. FigureQJa) toQJc) 
show the results of the system approaching the phase 
boundary from the SDW side, as the nearest-neighbor 
interaction V increases from 1.0 to 4.5 (while keeping 
fixed U = 10). For V = 1.0 and 3.0 the response of 
the system shows only a slight increase of mostly short- 
distance correlations during and after the pulse. In con- 
trast, at V = 4.5 the system displays clear tendency 
towards a sustainable charge order enhancement that 
emerges as the external pulse reaches its maximum and 
remains nearly constant after the pulse. By performing 
a time-dependent density-matrix renormalization group 
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FIG. 1: (color online) Time dependence of the charge-charge 
correlations as functions of distance (labeled by j) for 14-site 
lattice. The laser pulse with Gaussian magnitude modula- 
tion reaches its full strength at t = 12.5, as indicated by 
solid lines. Without exception, the pumping frequencies are 
set to match the resonance peaks of the optical absorption 
spectrum. Parameters: U = 10, St = 0.02, M = 100. (a) 
V = 1.0, w pump = 7.1, A = 0.10; (b) V = 3.0, w pump = 6.1, 
A = 0.30; (c) V = 4.5, o: pump = 4.0, A = 0.07; (d) V = 5.5, 
Wpump = 4.1, A = 0.60. 



method under open boundary conditions on larger sys- 
tem size, up to 30, similar behavior can be observed (not 
shown here). We note that the length scale with which 
prominent charge order enhancement can be identified in 
the case of V = 4.5 is confined within distance less than 
ten, in the units of lattice constant. The opposite effect 
is found when starting from the CDW side of the phase 
diagram, at V = 5.5 [Fig. HKd)], where after the pulse, 
the charge order is substantially diminished, along with 
a limited recovery of SDW order. 

For deeper understanding of the photoinduced charge 
order enhancement we note that the quantum phase tran- 
sition between SDW-CDW phases is driven by the com- 
petition between the energy cost for doublon genera- 
tion and the energy gain due to the attraction between 
doublon- holon pairs. The transition point is roughly 
V « £7/2. In the SDW phase, we expect charge-order 
favorite states to proliferate in the low energy regime as 
the transition point is approached, which makes them 
likely candidates to be picked up by the laser pulse. This 
is in contrast to the case when the system is far from the 
transition point where CDW states represent high-energy 
states. For a more quantitative analysis, we compare the 
spectra of systems with different values of V in the time- 
independent Hamiltonian, i.e., Eq. ([T]) but A(t) = 0. In 
order to describe the CDW order with finite spatial ex- 
tension, we define 

L-l L c 

Ocdw = 7^ E E(" 1 )' ~ !) ^ - !) ■ ( 5 ) 

C i=0 3=1 

Here, L c is introduced as a cut off parameter for the 
correlation length. The expectation values of the order 
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FIG. 2: (color online) The expectations of CDW order param- 
eter of eigenstates for 10-site lattice with V = 1.0, 3.0, 4.5, 
and 5.5. The energy E is measured from the ground state. 
Only the data of states with P to t = Stot = are shown, up to 
E = 60. Note that the largest value of (Ocdw) in (a) through 
(d) is marked by a cross. 




parameter (Ocdw) for eigenstates of a smaller 10-site 
system with L c = 5 are produced in Fig. [5] Comparison 
between different excited states shows that we can dis- 
tinguish states with large values of the CDW order. Note 
that despite the external laser pulse, the total momentum 
Ptot and the total spin S t ot remain good quantum values 
throughout the time evolution. For the ground state of 
a 10-site lattice at half filling, we have P tot = Stot = 0. 
Accordingly, only results of the eigenstates with the same 
values of Ptot and Stot are displayed. 

In Fig. [21 we find that in the case of small V = 1.0, 
the eigenstates with predominant CDW features are posi- 
tioned in the high energy regime, located around E ~ 45 
above the ground state energy. With the increase of V 
the states with large values of (Ocdw) move towards the 
lower part of the energy spectra. At V = 5.5 the eigen- 
state with (Ocdw) ~ 0.65 turns to be the ground state 
[Fig. [5[d)]. Keeping this picture in mind, it becomes 
more plausible that a well-tuned laser pulse may trigger 
the enhancement of the charge order on the SDW side of 
the phase diagram. Our numerical calculations suggest 
that the necessary precondition for such enhancement is 
the proximity of the system to the phase boundary as 
well as matching conserved quantum numbers between 
the SDW and CDW phases. 

We further elaborate on the condition for the emer- 
gence of the CDW order enhancement induced by the 
laser pulse. To this effect we perform parameter-sweeping 
calculations on the 10-site lattice. We sweep the pump- 
ing frequency w pump and laser intensity Aq and carry out 
the time evolutions up to t = 52.5. For a given pair of 
w pump and A , we then calculate the expectation of CDW 
order (OcDw)av, and the energy increase measured from 
the ground state AE, by averaging on the last 50 time 
steps (corresponding to time length At = 5). Contour 
plots of (OcDw)av and AE arc shown in Fig. [31 To facil- 



FIG. 3: (color online) Contour plots of final time-evolution 
results of CDW order (OcDw)av (first column) and energy 
increase AE (second column) as functions of u; pump and A , 
obtained by averaging on the last 50 time steps (time length 
At = 5), for 10-site lattice. Here, we take St = 0.1, M = 30. 
(a), (b) V = 1.0; (c), (d) V = 3.0; (e), (f) V = 4.5; and (g), 
(h) V = 5.5. 



itate the analysis we present in Fig. 2] the corresponding 
optical spectra, obtained from the imaginary part of the 
dynamical current-current correlation function: 



lmxj(u>) = jYl 



(n\j\0) 5(uj-E n + E ). 



(6) 



where |0) and \n) represent the ground state and excited 
states with energy Eq and E n , respectively. The current 

operator j defined as j = -i Y^%, a ( c !,o- c i+i,ff - H.c.^ . 

In Fig. [31 we notice the emergence of similar patterns 
between (OcDw)av (left) and AE (right) in each row 
panel. When starting from the SDW side [Figs. [UJa)- 
EUf)], in the case of moderate values of Aq, the contour 
plots representing the CDW enhancement and energy in- 
crease roughly match. On the other hand, starting from 
the CDW ground state, when V = 5.5 [Figs. [3jg) and 
[2[h)], the effect is just the opposite: the energy increase 
occurs along the destruction of charge order. Further- 
more, the positions of the enhanced energy stripes largely 
coincide with the resonance-peak positions of optical ab- 
sorption, as presented in Fig. [H except for some cases in 
the large Aq region, which will be discussed later. From 
these observations, we can conclude that at moderate 
values of Aq, in order to obtain enhanced (OcDw)av, the 
incoming photon frequency should be tuned close to ab- 
sorption windows, that match resonance peaks of the op- 
tical spectrum. 
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FIG. 4: (color online) The imaginary part of the dynami- 
cal current-current correlation function for 10-site lattice with 
V — 1.0, 3.0, 4.5, and 5.5, obtained by Lanczos method. The 
S function is approached by a Lorentzian with a width of 0.1. 



Let us investigate the case of V = 4.5 in more detail. 
The first optical peak is located at 10 ~ 4.3 [Fig. Sfc)]. 
States with enhanced CDW order are located around the 
same energy (the highest one with energy 3.9), as shown 
in Fig. [H[c). Not surprisingly, a sustainable enhance- 
ment of (OcDw)av can be found with w pump — 4.3 when 
Aq ~ 0.08 [Fig. EJe)] . On the other hand, systems with 
V = 1.0, and 3.0 [Figs. [3^ a) andEJc)] display no, or at 
most a very slight, increase of (OcDw)av around the res- 
onance frequencies, which match the optical frequency 
peaks, i.e., u ~ 7.5 and 6.3, respectively [Figs. 0|a) and 
0Jb)]. We thus propose two conditions for the observa- 
tion of the enhancement of (OcDw)av after the pulse: 1) 
the current matrix element (n|j|0) in Eq. (|6]) should be 
sizable and 2) (n\OcBw\n) should be large as well. 

Further detailed analysis on the V = 4.5 case shows 



that at Wpump ~ 4.3, with Aq growing beyond 0.1, AE 
temporally drops, reaches a local minimum and then 
keeps increasing [Fig. |3£f)]. The CDW signal up to 
Aq ~ 0.25 follows the same trend and then decreases 
[Fig. EX )] - This suggests a third condition to obtain en- 
hancement of (OcDw)a V : the pulse should deliver the 
optimal energy increase which can be controlled cither 
by Aq, or another parameter t^. In some simulations 
we have noticed a temporary increase of CDW order 
that lasts only up to the midpoint of the pulse action 
and then diminishes. The reason is that the system has 
passed the CDW enhancement region during the pro- 
cess - it has gained excessive energy Such a case can 
be found in Fig. QTb) (V = 3.0, Aq = 0.30), where we 
observe that two successive enhancements of the CDW 
signal appear during the pulse and then die out. The 
fact that the double-peak structure also emerges in the 
corresponding energy profile is consistent with the above 
argument. Moreover, during the irradiation, photons can 
be absorbed or emitted at different stages following the 
time evolution. 



We also notice the enhancements of the CDW signal 
with the accompanying energy increase when w pump is 
less than the resonance frequency, such as for V = 4.5, 
when Wpump « 2 [Fig. EJe)] . We speculate that this ef- 
fect is generated by the multiphoton, or more precisely, 
two-photon process, partly supported by the fact that it 
appears at large Aq and the energy increase is close to 
4, the same as what happens when the w pump is around 
4. The analysis on this possible multiphoton process in 
the off-resonant regime is out of reach of the present Let- 
ter 0. 

In conclusion, we address the issues of the thermo- 
dynamic limit and feasible experimental realizations of 
the proposed CDW enhancement. It is well known that 
for the ID extended Hubbard model, excitonic bound 
states can be found at the edge of the optical spectrum 



when V > 2t h 



The existence of well-formed cx- 



citons is a preliminary condition of the enhancement of 
the CDW order parameter after a laser pumping. With 
further increase of V, we can expect that at the edge 
of the absorption spectrum, charge-order states multiply 
and gradually dominate. They can thus be captured by 
a laser pulse with well tuned frequency and strength, as 
proposed by our numerical simulations. 

Although it is quite difficult to find a quasi- ID Mott 
insulator with proper U and V values that position the 
system near the phase boundary, there are other mecha- 
nisms which can take up the role of V. When additional 
on-site Holstcin phonons arc taken into account, it has 
been shown that the system can be driven from SDW 
to CDW phase by electron-phonon interactions [2(| [27j ■ 
Representative materials near the phase boundary with 
substantial electron-phonon interaction can be found in 
tetracyanoquinodimethane (TCNQ) series [28|. Another 
way to approach the phase boundary can be chemi- 
cal substitutions, as in halogen-bridged transition metal 
compounds [29| . Moreover, a transient phase transition 
from CDW to Mott or metallic phases has already been 
observed [3(j. These might prove to be useful in the 
search of proper candidates to observe the proposed ef- 
fect. As an alternative to electronic systems, CDW en- 
hancement can possibly be realized in cold atoms |3l| . 

Finally, we would like to emphasize that, while our 
work has concentrated on a particular CDW order en- 
hancement induced by a laser pulse, the proposed mech- 
anism may be effective in other cases when observing the 
response of correlated systems to an external perturba- 
tion near a phase boundary such as enhanced antifcrro- 
magnetic correlations, stripes, or possibly even supercon- 
ducting fluctuations. 
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